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We solve the neutron-deuteron (nd) scattering in the Faddeev formalism, employing 
the AW sector of the quark-model baryon-baryon interaction fss2. The energy dependence 
of the NN interaction, inherent to the resonating-group formulation of the quark-model 
£SJ ■ baryon-baryon interaction, is eliminated by the standard off-shell transformation utilizing 

the normalization kernel of two-cluster systems. An extra nonlocality originating from this 
procedure is very important to reproduce all the nd scattering observables below E n < 65 
MeV. A new algorithm to solve the Alt-Grassberger-Sandhas (AGS) equations is presented 
in the Noyes-Kowalski method. The treatment of the moving singularity from the three-body 
l/") ' Green function is shown in detail, which is based on the spline function method originally 

developed by the Bochum-Krakow group. The predicted elastic differential cross sections 
reproduce the observed deep cross section minima at E„ = 35-65 MeV and (9 cm = 130° -135°, 
which is consistent with the nearly correct triton binding energy predicted by fss2 without 
the three-body force. 

Subject Index: 205 

i £h t . §1. Introduction 

The QCD-inspired spin-flavor SUq quark model (QM) for the baryon-baryon 
^ . interaction, developed by the Kyoto-Niigata group, is a unified model describing 

If} | interactions between full octet-baryonsP^ 1 It is given by the Born kernel formulated in 

the resonating-group method (RGM) for interacting three-quark clusters. The short- 
range part is described by an effective one-gluon exchange, while the medium- and 
long-range parts are dominated by meson-exchange potentials between quarks. The 
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■ model parameters are constrained to reproduce all the two-nucleon data and available 

low-energy hyperon-nucleon scattering data. These QM baryon-baryon interactions 
are characterized by the nonlocality and energy dependence inherent to the RGM 
framework. For example, the Pauli-forbidden state appears on the quark level in 
certain channels of the strangeness sector, as a result of the exact antisymmetrization 
of six quarks. The short-range repulsion described by the nonlocality of the quark- 
exchange kernel gives quite different off-shell properties from the standard meson- 
exchange potentials. The energy dependence of the interaction is eliminated by 
the standard off-shell transformation, utilizing a square root of the normalization 
kernel N™'** This procedure yields an extra nonlocality, whose effect was examined 
in detail for the three-nucleon (3N) bound state and for the hypertritonP The 
advantage of the larger triton binding energy caused by the QM nucleon-nucleon 
(NN) interaction, namely, the deficiency of 350 keV, predicted using the most recent 
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model fss2, is still much smaller than the standard values of 0.5 - 1 MeVP given by 
the modern meson-exchange NN potentials. It is therefore interesting to examine 
predictions based on the QM NN interaction for the 3iV scattering, especially in 
this renormalized framework without the explicit energy dependence of the RGM 
kernel. 

An attempt to investigate the effect of nonlocality of the underlying two-body 
interaction on three-cluster systems composed of composite particles is actually not 
new. Even if we restrict our interest to the 3iV system, there are many investiga- 
tions to improve the triton binding energy and the neutron-deuteron (nd) scattering 
observables, from various types of nonlocal AW interactionP'~E3J There is, however, 
no extensive investigation of the 3iV system, employing a consistent NN nonlocal 
interaction that reproduces all the two-nucleon data with satisfactory accuracy. 

In this series of papers, we apply our QM NN interaction fss2 to the nd scat- 
tering in the Faddeev formalism for systems of composite particles EJ-O The 
Alt-Grassberger-Sandhas (AGS) equations-^ 1 are solved in the momentum repre- 
sentation, using the off-shell RGM t-matrix obtained from the energy-independent 
renormalized RGM kernel. The Gaussian nonlocal potential constructed from the 
fss2 is used in the isospin basis The singularity of the NN i-matrix from the 
deuteron pole is handled by the Noyes-Kowalski methodP^ 1 '^ To the best of our 
knowledge, this simple method to solve a singular Lippmann-Schwinger equation has 
never been seriously applied to the AGS equation, despite the fact that this method 
was originally developed for application to three-body problems. We will show a 
new practical algorithm to solve the AGS equations in this Noyes-Kowalski method. 
Another notorious moving singularity of the three-body Green function for the free 
motion is treated by the standard spline interpolation technique developed by the 
Bochum-Krakow groupP^~^2J Here again, we will give a detailed procedure and 
formulations since they do not seem to be available in the literature. We mainly 
use the channel-spin formalism, which is convenient to discuss the nd scattering. 
The accuracy of the numerical calculations can be checked by examining the optical 
theorem for the nd scattering. The NN interaction up to / max = 4 with a suffi- 
cient number of partial waves of the three-body system is included, resulting in the 
maximum neutron incident energy of about E n = 65 MeV in the laboratory system. 
Here, / max is the maximum value of the two-nucleon angular momentum included 
in the calculation. In this paper, we only show the results of total and differential 
cross sections of the nd elastic scattering. Polarizations and deuteron breakup cross 
sections are discussed in subsequent papers@ The details of the Gaussian nonlocal 
potentials and an application to the S-wave nd scattering length will be given in 
separate papers.® 

The organization of this paper is as follows. In §2.1, we first recapitulate our 
QM baryon-baryon interaction and then outline the whole procedure to solve AGS 
equations for the nd scattering. An emphasis is put on how to calculate an extra 
nonlocal kernel originating from the elimination of the energy dependence in the 
RGM kernel. The essential part of this analytic formulation is given in Appendix A. 



A preliminary report on the present subject is found in Refs. l23p ~t25 p . 
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In §2.2, we exhibit a new formulation of AGS equations, incorporating the deuteron 
singularity of the NN f-matrix in the Noyes-Kowalski method. The basic equation 
is convenient to derive the optical theorem by a simple manipulation in the operator 
formalism. The contribution of the breakup cross sections to the optical theorem 
is thoroughly discussed in §2.3. Practical calculations are made by the procedure 
given in §2.4, where the moving singularity of the three-body Green function for 
the free motion is processed by the subtraction method and the spline interpolation 
technique. The detailed procedure to calculate the basic integral Qk^v is given in 
Appendix B. The total cross sections derived from the optical theorem and nd elastic 
differential cross sections are shown in §3.1 and §3.2, respectively. It is found that 
the predicted elastic differential cross sections reproduce the observed deep cross 
section minima on the high-energy side, which is consistent with the nearly correct 
triton binding energy predicted using fss2 without the three-body force. The last 
section is devoted to a summary. 

§2. Formulation 

2.1. Quark-model baryon-baryon interaction 

The QM baryon-baryon interaction is a low-energy effective model that intro- 
duces some of the essential features of QCD (quantum chromodynamics) character- 
istics. The color degree of freedom of quarks is explicitly incorporated into the non- 
relativistic spin-flavor SUq quark model, and the full antisymmetrization of quarks 
is carried out in the RGM formalism. The gluon exchange effect is represented in 
the form of the quark-quark interaction. The confinement potential is a phenomeno- 
logical r 2 -type potential, which has a favorable feature in that it does not contribute 
to the baryon-baryon interactions. We use a color analogue of the Fermi-Breit (FB) 
interaction, motivated from the dominant one-gluon exchange process in the high- 
momentum region. We postulate that the short-range part of the baryon-baryon 
interaction is well described by the quark degree of freedom. This includes the 
short-range repulsion and spin-orbit force, both of which are successfully described 
by the FB interaction. On the other hand, the medium-range attraction and long- 
range tensor force, especially afforded by the pions, are extremely nonperturbative. 
These are therefore most relevantly described by the effective meson exchange po- 
tentials (EMEPs). The most recent model fss2 includes the vector- meson exchange 
EMEP, in addition to the scalar- and pseudoscalar-meson exchange potentialsP^ 1 '^ 

The quark-model Hamiltonian h consists of the phenomenological confinement 
potential U-j , the colored version of the full FB interaction Uf B with explicit quark- 
mass dependence, and the EMEP UfJ generated from the scalar (i?=S), pseudoscalar 
(PS), and vector (V) meson exchange potentials acting between quarks: 
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The RGM equation for the relative-motion wave function x( r ) reads 

(<j>(3q)<t>(3q)\E - h\A{<l>(3q)<j>(3q)x(r)}) = , (2-2) 

where 4>(3q) is a simple harmonic-oscillator (h.o.) shell-model wave function for the 
three-quark clusters. We solve this RGM equation in the momentum representa- 
tion.^ If we rewrite the RGM equation in the form of the Schrodinger-type equa- 
tion as [e — ho — VRGM(^)]x( r ) = 0, the potential term, Vrgm(£) =Vd + G + eK, 
becomes nonlocal and energy-dependent. Here, Vb represents the direct potential of 
EMEPs, G includes all the exchange kernels for the interaction and kinetic-energy 
terms, K = 1 — N is the exchange normalization kernel, and e is the total energy in 
the center-of-mass (cm) system, measured from the two-cluster threshold. We cal- 
culate the plane- wave matrix elements of Vrgm(£), and set up with the Lippmann- 
Schwinger equation for the RGM i-matrix. This approach is convenient to proceed 
to the G-matrix calculations^ 1 >E3 anc i i Q the Faddeev calculations with some special 
considerations of the Pauli-forbidden statesP^ 1- ^ 

The energy-independent renormalized RGM kernel y RGM for a two-cluster sys- 
tem is given byS 

yRGM = y D + q + w _ ^.3) 

The nonlocal kernel W appears through the elimination of the energy dependence, 
and is given by 

W = A^=h^=A - h . (2-4) 



Here, N = 1 — K is the normalization kernel, h = ho + Vb + G with ho being 
the kinetic energy for the two-cluster relative motion, and A = 1 — \u)(u\ is the 
two-cluster Pauli projection operator, where \u) is a Pauli-forbidden state satisfying 
K\u) = \u). The Born kernel W(p,p') of Eq. (|2-4p . or its partial wave component, 
is calculated accurately using an analytical procedure to obtain JC = A{\/\fN — \)A 
in the momentum representation. This is discussed in Appendix A. In the NN 
sector, no Pauli-forbidden state appears on the quark level, so that we can simply 
set A = 1 in the following formulations. An advantage of using 

yRGM ig that the twQ _ 

cluster RGM equation takes the form of the usual Schrddinger equation in the Pauli- 
allowed model space, and the relative wave function is properly normalized. This 
Schrodinger-type equation for the relative wave function gives the same asymptotic 
behavior as the original RGM equation, thus preserving the phase shifts and physical 
observables for the two-cluster systems. The difference between the previous energy- 
dependent RGM kernel, V KGM {e) = Vq+G+eK, and V RGM in Eq. is essentially 
a replacement of A{eK)A with W . The value of e is, however, not properly defined in 
the three-cluster system, in particular, for the scattering systems. In the following, 
we will consistently use the energy-independent renormalized RGM kernel V M in 
Eq. (|2-3p . both for the bound-state solution and the scattering problems. 
The three-cluster equation for the 3iV bound state is written as 

[ E - Ho - V* GM - V> RGM - F 7 RGM ] $ = o , (2-5) 
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where a, f3, and 7 denote three independent pairs of two-cluster subsystems, Hq is 
the free three-body kinetic-energy operator, and V^ GM stands for the RGM kernel 
in Eq. (|2-3|) for the a-pair, etc. In Ref. EJ), we have solved Eq. (|2-5|) in the Faddeev 
formalism 

i/j = G tPip , (2-6) 

where Go = Gq{z) = l/(z — Hq) is the three-body Green function for the free 
motion, P = P(i2)P(l3) +-f(i3)-f(i2) * s a permutation operator for the rearrangement, 
and t = t{z — ho) is the AW i-matrix derived by solving the Lippmann-Schwinger 
equation t = v + vG$t with v = y RGM in Eq. (|2-3|) . Here, z = E is the total energy, 
which is below the deuteron energy (< 0) for the 3iV bound state. The energy 
argument of t is z — ho, where ho is the kinetic-energy operator for the relative 
momentum q = (2k-^ — ki — k,2)/'i. Another momentum for the two-nucleon relative 
motion is denoted by p = (k\ — k2)/2 with ho being the corresponding kinetic-energy 
operator. The three-body kinetic-energy operator is therefore Ho = ho + ho- I n Eq. 
(|2-6|) . ijj is the Faddeev component that yields the total wave function ^ in Eq. (|2-5[) 
through & = (l + P)ip. The basic equation for the nd scattering is the AGS equation, 
which can be expressed as^D 

U\</>) =G ~ 1 P\cf ) } + PtGoU\ ( i ) } , (2-7) 

where \4>) = \qo-,^d) is the plane- wave channel wave function with |^) being the 
deuteron wave function. In this case, z = E + iO is the total energy approached from 
the upper side of the real axis in the complex energy plane, and E = E cm + with 
-E^cm = (3h 2 /AM)qo 2 being the neutron incident energy in the cm system. We use an 
average nucleon mass M = (M p + M n )/2 in the isospin formalism. The scattering 
amplitude for the elastic scattering is obtained from (<p\U\(p). It is important that 
Eq. (|2-7p also provides information on the full breakup process of the deuteron. The 
transition amplitude for the breakup process is given by 

C/ o |0) = (1 + P)tGoU\<j)) = (1 + P)T\<j)) , (2-8) 

where T = tGoU corresponds to the three-body f-matrix. The breakup cross sections 
are obtained from the amplitude (pg|£/o|<A) with the corresponding energy E = 
(h 2 /M)(p 2 + (3/A)q 2 ). 

2.2. Noyes-Kowalski method for the singularity of the NN t-matrix 

For the description of the nd elastic scattering, it is convenient to use the channel- 
spin representation, in which the angular-spin functions are defined through 

|p,q;123) = 7) <7|P,Q;123> 

7 

<p,?;123|7> = [[y A (p) Xs< (l,2)] 7 xi i(3)] J „ , , (2-9) 



' 2 ' 



with 7 = £[(\s)I±]S c ; JJ Z ; {t±)\T z . Here, Xs*(l,2) etc. are the spin-isospin wave 
functions. In particular, the deuteron channels are specified by 

ld = 4(Al)li]S c ; JJ Z ; (01)17; , (2-10) 
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with A = and 2, corresponding to the S- wave and D-wave components, respectively. 
We use t = |7}ty(7|, and separate the 7-sum into 7 ^ (7^) and 7 G (7^). The 
i-matrix of the deuteron channel has the deuteron pole, which we explicitly separate 
as 

t ld = — *a . (2-n) 

Z — <lQ ld — £d 

If we use the completeness relationship, J °° q 2 dq \q){q\ = 1, and the spectral de- 
composition of the two-nucleon Green function, we can easily show 



z - h Qld 

where we have used a notation 

q M 



AM f°° t \n\ ^ 

- = ^ V q 2 dq\q)-^^(q\-ic(\ q o)t ld [ qo ](qo\) , (2-12) 
£d 6n Jo Qo — Q 



2vr- 



3h 2 ' 



hM = hd (4^0 ~q) + £d 

haiqo] = t-y d ( £ d) = goiedT^di^chgoiedT 1 , (2-13) 

with go(sd)~ 1 \'4 , d} = {^d — ho)\ipd) = v\tpd) being the deuteron solution. For the 
separable deuteron residue, we in fact need to make a distinction between the S- 
wave and D-w&ve channels by writing t 7d (e^) as 

t-y d ,iA £ d) = go( £ dy 1 \iP ld )(iP%\go(edr 1 , 

m^d) = E<v>7> 7d > = 1 , £ hd)\^ d ) = Hd) . (2-i4) 

id id 

In Eq. (|2-12p . we should note 

\qo)h d { £ d){m\ = go{edT 1 \qo,tPd){qo^d\go{£dT 1 = G , o _1 |^)<0|G'o~ 1 , (2-15) 

since 

go{ £ d)~ l \qo,Tp~ ld ) = (£d - h )\qo,ipd) = ( 4^9o 2 + e d - h - h J \q ,ip d ) 

= Go^lqo^d) = Go- 1 ^) , (2-16) 
with \4>) = \qo,ipci)- In fact, a more rigorous expression is 

hd)\qo)t-y d ^ d (£d)(qo\(ld\ = £ l7d)G'o _1 ko,V'7d)(9o,^7JGo~ 1 (7d| 

ld,ld 7d,7d 

= Go-^o^^^c^lGo" 1 = Go-VX^IGo -1 • (2-17) 



Quark-Model Baryon-Baryon Interaction Applied to nd Scattering 



7 



Here, ■jd is defined from 7^ in Eq. (|2T0p by just replacing A with A. Since this tensor 
coupling is cumbersome in the complicated expressions, we shall use in the following 
a simple notation to represent t~ d ~ d by i 7d . After all, we have obtained 



^ *- ho ld - e d 3H 2 ^ J q 2 - q 2 

-icGo-^i^Go- 1 . (2-18) 



We write the first term of the right-hand side (r.h.s.) as t 7d and generalize the 
expression to all the channels: 

t, = {. f I , fa lihd) . (2.i 9) 

This yields a separation of the two-nucleon singularity from the full i-matrix: 

t = t-icG o ~ 1 \0)(0\G o - 1 , (2-20) 
where t = £ 7 |7)*7<7l (= £ 7 ^( 7d ) MMtI + E 7d I 7rf)*7d(7rf|) with 

= S V r^dq \q,l d )-^^(q,7d\ • (2-21) 

•J" jo Qo — q 

The difference between i and t appears only when q = qQ in the deuteron channel, 
and t also satisfies the same basic t-matrix equation t = v + vG^t. If we use Eq. 
(I2 : 20l) in Eq. ([22]), it becomes 

U\<f>) = Go^P^) [1 - ic{<f>\U\<f>)] + PtG o C/|0) , (2-22) 

since PG _1 = Go -1 P. If we define U'\<f>) by 

I7|0) = l7»[l-ic(0|tf|0>] , (2-23) 

it satisfies the following equation: 

V\<f>) = G Q ~ X P\4>) + PtG Q U'\<j>) . (2-24) 

We set Z = {4)\G Q - l P\(f)) and multiply Eq. ||2^3|t by P\</>) Z' 1 (<j)\ from the left-hand 
side (l.h.s.), and obtain 

P\4>) Z- 1 {<W\4>) = P\4>) + P\4>) Z- 1 {<j>\PtG Q U'\ct>) . (2-25) 

We further multiply Eq. (|2-24p by Go from the l.h.s. and subtract Eq. (|2-25p from 
the result. Then, the first terms of the r.h.s. cancel, and we obtain 

G Q U'\<j>) = P\<f>) Z- 1 (<f>\U'\<f>) + WtG o U'\0) , (2-26) 

where we have defined W by 

W = G P-P\(j)} Z- 1 {<j>\P . (2-27) 
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Thus, if we set 

G Q U'\<t>) = Q\4>) Z- 1 (4>\U'\4>) , (2-28) 

Q\4>) satisfies 

Q\4>) = P\4>) + WtQ\4>) . (2-29) 

We should note that W in Eq. (|2^2Tj) satisfies {(j)\G ~ l W = and WGo _1 |0> = 0. 
Thus, if we multiply Eq. (I2 : 29ll by (0|G O _1 from the l.h.s., we obtain 

(0|Go _1 O|0> = (^IGo-^l^) = Z , (2-30) 

which is consistent with the definition of Q\(f>) in Eq. (|2-28|) . Furthermore, i 7 in Eq. 
(|2-29p can be restored to i~ for the deuteron channel 7 = 7d owing to WGo' 1 ^) = 0, 
thus allowing us to write t as t. 

We should note that the channel wave function \<f>) is actually \4>;{IS C )J J z ) 
in the channel-spin representation with at most three possible configurations, i.e., 
(£S C )J = (J + f , I) J, (J - |, |)J, and (J - i, a) J for the parity vr = (-l) J ~K 
and (J - |, I) J, (J + I, ±)J, and (J + ±, |)J for vr = (-l) J+ i Namely, we should 
understand 

\(f>;(£S c )JJ z ) = \q Q ,tp d ;(iS c )JJ z ) = P 2 d P \P,Qo,ld) ip-y d (p) > ( 2-31 ) 

a Jo 

using 7^ defined in Eq. (|2-1U|) . The basic relationship in Eq. (|2-3Up implies that Q\4>) 
is a modification of P\4>) by the effect of the nonsingular interaction t, and the real 
symmetric matrix Z in Eq. f|2-30j) plays an essential role in the following discussion. 
In particular, Z is a nonsingular matrix as seen in Eq. 42381)0 Thus, if we define 

Q\4>) = Q\<j>) Z- 1 , %) = P\<j>) Z- 1 , (2-32) 

and multiply Eq. (|2-29j) by from the r.h.s., we obtain our final equation 

Q\<f>) = P\4>) + WtQ\cf>) . (2-33) 

In order to derive the scattering amplitude, we multiply Eq. (|2-24p by (cp\ from 
the l.h.s. and use Eq. (|2-28p . Then, we find 

(<P\U'\<P) = Z + (0|PtQ|0) Z- 1 {(p\U'\4>) . (2-34) 

Then, if we set 

= Z- 1 (<j>\PtQ\4>) Z~ l = (<j>\PiQ\4>) , (2-35) 

we find 

{<f>\U'\<f>) =Z [l + (0|X|^)(0|C/»] = [Z- 1 - WXW]- 1 . (2-36) 

*' This can be proved if g ld .-y' d (<?o, qo \ x) in Eq. (|2-41jl is a positive- or negative-definite function, 
which is true at least for the dominant 5-wave deuteron component. 
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This expression, together with Eq. (|2-23|) . yields 



^\U\4>) = [l + ic{<j>\U' 



-1 



\u' 



\X\6)+ic-l\ 



(2-37) 



A prescription of the principal- value integral still remains in Eq. (|2-35[) for the 7 = 7d 
term. The residue at q = qo is 



^2 1o 2 (4>\P\qo, ld)% d [Qo] (Qo, ld\Q\4>) 



id 



^2 qo 2 (<f>\P\qo,jd)go(£d) 1 \^- yd }(^y d \go(£d) l (qo,j d \Q\ ( l ) ) 

q 2 ((ftPGo-^MGo^Qlt) = q 2 Z 2 , (2-38) 



where Eq. (|2-30p is used in the last part. Thus, we only need to subtract Eq. (|2-38p 
since f Q °° dq l/(q 2 - q 2 ) = 0: 

(ct>\x\<i>) = (^\ p h) h <7lW> 



+ w> dq 



go 2 - q 2 



q 2 ^2(<p\P\q,7d) ty d [q] {q, Jd\Q\<j>) ~ 9o 2 



7rf 



(2-39) 



We should note that the final equation (|2-33p still contains a delta function from 
the rearrangement factor 



{p,q,i\P\p',q',7' 



dx 



5(p-pi) 



p 



X+2 



g 1;1 >(q,q ,x) 



, J{p' -P2) 



P 



,'X'+2 



(2-40) 



Here, p\ = p(q' , q/2; x) and P2 = p(q, q' /2; x) with p(a, b; x) = Va 2 + b 2 + 2abx. We 
use the spline interpolation for the variable p and introduce the weight factors for 
the Gauss-Legendre quadrature. For these, we use the notations 



P^uk = p{q^ qu/2, x K ) = q/j 2 + q v 2 /A + q~^x~ K , 

Si^K = Si {Pijuk) , (coefficients for spline interpolation) 



/ 1 \ A 1 

Bjfj.,j'u,K = qii\f^\~iqu\f^u ( ) 7; 



2^7,7' qu\ %k) 



PfJUK 



Uk , (2-41) 



where q^, x K , etc. are Gauss-Legendre discretization points and u) K , etc. are 
the corresponding weights. The two-nucleon relative angular momenta, A and A', in 
BjUjj'vjK are related to 7 and 7', respectively. For the on-shell momentum q^ — > qo, 
we assume qoy/uio = 1. Similarly, for k = 0, appearing in the pole prescription for the 
x-integral later, we assume ujq = 1. The weight factors are introduced, for example, 
as 



W, 



qfjju~vP 2 Ui {puqfj^\W\pj,q u ,i) p 2 Uj q v ^JZj, 



(2-42) 
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The key relationship to avoid the appearance of the delta function for p is the 
replacement 

ujj 5 (pj - p^ UK ) = S jflVK , (2-43) 
which can be proved from the integration formula 

r 6{p/ np') = f( P , UK ) = y: s jlWK fto) , (2-44) 

Jo P j 

for arbitrary smooth functions f(p). For the spatial part of the deuteron wave 
functions, F\(p) = ip~ fd (p), we also use the notation F,j\ = F\(p) = F\(pj)/pj , but 
apply the spline interpolation not to F\(p) directly, but to (p\go~ l {£d)' l Pd) ■ This is to 
avoid the numerical inaccuracy of the spline interpolation and guarantee the exact 
relationship WGo~ 1 \4>) = P\4>) — P\4>) = 0. Namely, we define 

9jX = (Pj\9o~\£d)^d) = (e* ~ ^pfj F jX , (2-45) 



and calculate F\^ K through 



F\^ K = hd - JjPIok) ^5*iM0« 9j\ , (2-46) 

3 

instead of F\^ K = J2j s j^0n FjX- Similarly, we calculate F XOuK through 

/ h 2 \~ 1 

Fxou. = hd- mPIvk) Yl S i°™ 3jX ■ (2-47) 

j 

Then, by using the definition, P im = ppuii q^y/U^l (pi,q^,-f\P\<j)), we calculate P^ 
and Z through 

Pin~i = ^ ] ^ ] SiQ^ K B^^ ^Q K F\i , 

z = /^QixPiOjd = Z^2[ £d ~ m P ° 0k ) ^ ^ A00k B ^o,y d o,K Fx>ook ■ (2-48) 

iX n ' A, A' 

For Pi^ = Pi^y Z~ l , we trivially obtain 

Up = 3ixPi0 ld = Z Z- 1 = 1 . (2-49) 

iX 

For the two-body i-matrix, we define 

^2 



€=w;l w s,>(^» 2 -* 2 >+^)n 



4tt 
(2~ 



^3 tlluxs) («• W S ((? ° 2 - q " ] + £d ) ' (2 ' 50) 
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which is simply denoted by i 7 



The discretization points for are divided into at least two or three regions, 
depending on the incident energy. For qo < with = y (AM/3h 2 )\ed\ (the 
neutron incident energy E cm = (3ft 2 /4M)qo 2 < \^d\), we only have the nd elastic 
scattering, and the behavior of the two-body i-matrix is 



t : real /* oo \ 



(2-51) 



For go > (E cra > \£d\), the three-body breakup is possible and the three-body 
Green function Gq has a pole. The behavior of the two-body f-matrix in this case is 



Qii ■ [0, \J qo 2 ~ K d 2 } [\/qo 2 ~ K d 2 ,qo] Qo [<?o,oo) 
t : complex real /* oo \ 



(2-52) 



so that qu = y qo 2 — «<? is the threshold momentum for the deuteron breakup. 
The three-body Green function pole is treated in accordance with the subtraction 
method discussed in Glockle et al.'s review paper 

With these preparations, we can now write down the expressions for the numer- 
ical calculations. The most important matrix elements are 



A 



( T 



M 



x < 



log 



FO^I 



ft-q^qv 
for \x^ v \ > 1 , 

for \x 0fJ , v \ < 1 , 

(2-53) 



with xq^ = ((3/4)g M 2 - q 2 ~ qu 2 ) /(q^Qu)- The matrix elements ofW in Eq. fj2-42f) 
are calculated from 



Wi 



We can prove that 



, Yl W W7J0y' d 9j\> = 



We also define the main kernel for the linear equation by 



(2-54) 
(2-55) 

(2-56) 



For the numerical check, we extend the definition of M, 



M, 



5> 



in Eq. f)2-56f) to include 
(2-57) 
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for the on-shell 7' = 7^ and v = 0, where the residue of the off-shell i-matrix t^j in 
the deuteron channel is given by (a separable kernel) 



The basic AGS equation Eq. (|2-33p is then expressed as 



Pi, 



jwy' 



For Qi07 d calculated from 



(2-58) 



(2-59) 



vy i 



3 V 1 



we can show that 



Uq — 9i\ QiOld — 1 



(2-60) 



(2-61) 



similarly to Eq. (|2-49p . The explicit expression of Eq. ()2-39j) is 

Xis c 2'S> = {<t>£S c \X\(j)irs>) 



JM7 



4M 1 



^! ^M7d ^> Qimd 90 ^ 



7d,M 



(2-62) 



and the 5-matrix, SfL, s ,\ ™ Sc s = 5(.,i>5s a ,s' c — *2c U/g, s n ng\, is simply obtained by 
solving an equation 

[C^ -1 )^,.,^ ~~ Xis c #s' c + Sej'Ss^s',. Uy'sQW'S'j) = $W'fis e ,s'j ■ (2-63) 

2.3. Optical theorem 

The original AGS equation contains the full information of the unitarity for the 
three-body scattering. Here, we derive the optical theorem, starting from the AGS 
equation Eq. (j2-7[) or 



U = G - 1 P + PtG U , 



(2-64) 



for the systems composed of three identical particles. We first assume that there is 
no deuteron pole in the two-body i-matrix t. We take the hermitian conjugate of 
Eq. (p 7 ^ and replace G -1 P with U - PtG U: 

C/t = G _1 P + tfGjtfP = Go^P + E7+GoV6?o • Go -1 -? 

= Gq^P + U^GqU^GqU - rfGoWGoPtGoU . (2-65) 
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We again take the hermitian conjugate and use the fact that P and Go are exchange- 
able, 

U = Gq~ x P + U^GqUGqU - U j G j t^Go j PtG U . (2-66) 
The subtraction of Eq. (|2-66p from Eq. (|2-65p yields 

Z7+ - U = Z7 + G t (* t - t)G U - U^G q H\Gq - G ^)PtG U . (2-67) 
Here, we use a basic relationship for the two-body i-matrix, 

*t - t = tt(G t - Go)* (2-68) 

and 

G f - G = 2m5{E - H ) . (2-69) 

Then, we find 

U ] -U = 2-Ki U ] G H ] 5(E - flo)(l + P)tG U . (2-70) 

If we have a deuteron pole (or some equivalent divergence of the two-body t- 
matrix), we have to separate this divergence term using Eq. f|2- 20[) . Since t does not 
involve the divergence, we can apply the above discussion to U' in Eq. (|2-24p and t, 
and obtain 

C/'t -U' = 2-Ki U' ] G^5{E - H )(l + P)tG U' . (2-71) 

If we multiply Eq. (l2TT]l bv [l+ic(6\U\6)*](6\ from the l.h.s. and by \4>)[l-ic{(j)\U\(j))] 
from the r.h.s., the relationship in Eq. fj2-23|) yields 

(<j>\u\<j>)* - {<t,\u\<t>) = 2ic{4>\u\4>)*{<j>\u\<j>) 

+2m (<j)\U ] G Q W 5{E - H )(l + P)tG U\4>) . (2-72) 

We can derive another expression similar to Eq. (|2-72p . starting from 

[Z-'-^X^+ic] (d>\U\d>) = 1 , (2-73) 

in Eq. ([237|). We multiply Eq. (|275jl by (4>\U\4>}* from the l.h.s.: 

{^UW = mW [Z- 1 - (</>\X\4>) + ic] (<t>\U\<j>) . (2-74) 

We further take the hermitian conjugate of Eq. (|2-74p and subtract it from Eq. f|2- 74[) . 
Then, we find 

{4>\uw - {4>\u\4>) = 2ic (4>\u\4>)*(<t>\u\<t>) 

+{<t>WW MXW - (<P\x\<P)] (<W\<I>) ■ (2-75) 

Thus, the second term of Eq. (|2-75p corresponds to the second term of Eq. (|2-72p . 
This implies that the imaginary part of (<p\X\(j)) is related to the breakup cross 
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sections. We can also derive this from the basic equation in Eq. ()2-33p . We make 
{4>\P from Eq. (|2-33p and use this in Eq. (|2-35p . Then, we find 

(<I>\X\4>) = (<P\Q j tQ\0) - (^QWW^tQ^) . (2-76) 

We further take the hermitian conjugate of Eq. (|2-76p and subtract Eq. (|2-76|) from 
the result: 

{<t>\XW - {<t>\X\<t>) = (0|Qt(tt _ t)Q\4>) - (<f>\&V(W - Wl)tQ\<t>) . (2-77) 

Then, by using Eq. ([glgD for land - W = (G j - G )P derived from Eq. flMZD , 
we eventually obtain 

(0|X|0>* - (<j>\X\<j>) = (0|Q t ?(G o t - G )(l + P)tQ\^ 

= 2m (<f>\Q^5{E - H )(l + P)tQ\<j>) . (2-78) 

Note that the first term of Eq. (|2-77p gives the direct term of the breakup process, 
while the second term gives the exchange term. The last expression of Eq. (|2-78p 
corresponds to the second term of Eq. (|2-72p through Eq. (|2-75|) , since we have 

G U\<f>) = Q\4>)Z~ l {4>\U\4>) = Q\4>){4>W\4>) = Q\<P) • (2-79) 
After all, we have obtained the relationship 

(-2i) Im^U^) = 2ic J dq \{q Q , ^U^ 2 

+2m {<j)\&V5(E - H )(l + P)tQ\<f>) , (2-80) 

where the angular integral over q in the intermediate states is explicitly written. 

We should note that t in Eq. (|2-80p can be safely replaced with t, owing to the 
existence of the energy-conserving delta function of 5(E — Hq). In fact, if we use 
Eqs. dZSSD and (|2^30]l . we find 

tQ\<f>) =tQ\<t>)-icG - 1 \<t>)Z . (2-81) 

Here, the second term does not contribute to the last term of Eq. (|2-80p . since 
5(E - Ho)G - 1 \^) = 0. Thus, tQ\4>) in Eq. fl2~80]) can be replaced by tQ\<f>) = 
tGoU\4>) = T\(f>), where T = tGoU is the three-body t-matrix. If we further use 
(1 + pf = 3(1 + P) and the three-body breakup amplitude, U \4>) = (1 + P)T\4>), 
Eq. (|2-80p can be equivalently expressed as 

(-2i) Im(4>\U\4>) = 2^c J dq \(q ^ d \U\4>)\ 2 

+ 2m^JdpJdq5^E-^p 2 -^q^ \(p, q|[/ o |0)| 2 ■ (2-82) 

In the channel-spin representation, we take the spin sum for the initial spin states 
\S C S CZ ) and obtain 

(-2*) Yl iMfrS c s cz \u\frs c s cz ) 
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2ic jdq £ Yl \(qo,i>d;S' c ,S' cz \U\cP;S c S cz )f 



+2vui I dpi dq 5 (e - ^p 2 - ^qA £ £ |(p, g; S^IW; S C S CZ >| 



(2-83) 



We can carry out the partial-wave decomposition using 
{4>q f ;S' c S' cz \U\(j)q i ;S c S cz )= U(l>SQ,{JlSc) 

I'UJz 

xJ2(t'rn'S' c S'JJJ z )Y e , m ,(q f ) ^(ImSJS^JJ^Y^i) , (2-84) 
m' m 

with |qry| = |qfj| = go- Since the operator f7 is ./^-independent, we can prove 

£ (0; S c £* W; 5 C< S CZ ) = i- ^ (2 J + l)C/ ( J £5c)i( , 5c) . (2-85) 

In the first term of the r.h.s. in Eq. ()2-83j) . the q integral allows us to take the sum 
over S cz and J z in a similar way. In the second term of the r.h.s. in Eq. (|2-83|) . we 
first integrate over p and q, and change the S' c , S' cz sum to the 7 sum: 

1 f°° 2 , f°° 2 , r & 2 3h 2 

2m - p 2 dp q 2 dq5[E- —p 2 - —q 



Jo 



M AM 

Yl ^qv S cScz\U ] \p, q, j)(p, q, y\U \<f> Qi \ S c S cz ) . (2-86) 



X 

Sc->S C z T 



Here, the p-integral can be carried out from the 5 function. We use 



with po = p (q) = \J (3/4) (q M 2 - q 2 ) L which implies that q < q M = V 'qo 2 ~ ^d 2 - 
(This corresponds to the change from t to t before.) Thus, we find 

f°° ( h 2 3h 2 

p 2 dp q 2 dq 5 I E — —p 2 — tvtQ 2 
L \ M AMI 



/qm r°° f 2h 2 \ 

q 2 dq p 2 dp 5 {-JfPo(Po ~ p)j 



2h 2 

Thus, Eq. ()2-86p becomes 



M f qM , 2 

7^? dqq po- (2- 



1 M f gM _ ^ 

27Ti 3 2h? / dq q2p ° ^ ^ ^v ScSc ^ Uo \ po ' q '^ 

^(po,q,l\Uo\4> qi ;S c S cz ) . (2-89) 
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The rest is the same as in deriving Eq. (|2-85p . After all, the partial- wave decompo- 
sition of Eq. (|2-83p is given by 



E ( 2J + !) ImC/ (k),(^) = i^r E ( 2J + !) H>s> c ), 

(£S c )(e'SQJ 

/ ^<7 2 R) V (2J + 1)V | (po, q, j\U \q ,^ d ;(£S c ) JJ Z ) | 2 . 



(2-90) 

We multiply Eq. ff^QO]) by the overall factor (l/i)(3h 2 /q M), 
1 3h 2 

(- 2 )^^m E( 2J + 1 ) ImC/ (W c) = E (2^+i)i^),(^)i 2 

(« C )J {iS c ){t'S' c )J 



1 7T f qM 

+— 

47T 



/ dqq 2 p (2J+l)^|(po,(?,7l^oko,^;(^ c )JJ z )| 2 , (2-91) 

and use the elastic and breakup scattering amplitudes with a common factor 

f j 1 4M(2vr) 3 . _ *t4M j 

J(JPS' e ),(fi3 e )- (4^) 3n 2 (4tt) U (*'S'cU*Sc) - 2 3^" V'S&VSo) ' 

/ 7 W«) = -2^2 <Po,9,7l^o|ao,^; (i5 c )JJ«) • (2-92) 
Multiplying Eq. (|2 T 9T1) by the overall factor 47r{-(7r/2)(4M/3/i 2 )} 2 , we find 

^ £ (2 J + 1) Im/{ Sc) , (<Sc) = 4vr £ (2J + 1) |/ ( i 5 y )( , 5e) | 2 
(AS^J (£Sc){£' S' C )J 

+£ d? e ( 2J + 1) E i4S>)i 2 • ( 2 - 93 ) 

We find that the following relationship holds for each (iS c )J component: 

^ Im /(k),(^ c) = E \f(t>s' c ws c )\ 2 + i / ^ 9% E i/£)(9)P • 

(2-94) 

The relationship of the total cross sections is obtained by taking a spin average over 
the initial spin multiplicities ((2/ + l)(2s + 1) = 6 with 1 = 1 and s = |), namely, by 
taking the sum (4tt/6) X}(£s c )j( 2 ^ + 1) over Eq- (|2-94|) . For the practical calculation 
of the total breakup cross sections in Eq. (|2-94|) . it is most convenient to use the 
imaginary part of {(p\X\(p) through Eq. (|2-78p . It is rather straightforward to derive 



— lm f(ts c )MSc) ~ E \f(trs' ),{ts c )\ 2 

(i'SQ 

+ \n) M4^ f(£'S' c ),(iS c ) lmX (£'SQ,(e"S'J) rfff'S'MtSJ 1 ( 2 ' 95 ) 

v 7 w (i'S' c ),(e"s^) 
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where X^, s ,^ ^„ s ,^ is given in Eq. (|2-62p . 

2.4. Moving singularity of the three-body Green function 

The direct solution of Eq. (|2-59|) in terms of Eq. (|2-53|) causes a serious numerical 
problem, when it is applied to the energies above the three-body breakup threshold, 
namely, E cm = (3h 2 / AM)qo 2 > |e^|. If we use a restricted number of discretization 
points, say, 5 points for each interval of Eq. (|2-52|K it gives reasonable results. How- 
ever, if we increase the number of discretization points, the phase shift results of 4 S3 

2 

states, for example, deviate by more than several degrees. The origin of this patho- 
logical situation is the discretization points close to the logarithmic singularity in 
the kernel Eq. (|2-53|) . This is a notorious problem of moving singularities, appearing 
in any type of three-body model, that takes account of breakup processes. We avoid 
this by using the prescription given by the Bochum-Krakow group namely, 
applying the spline interpolation even to the q degree of freedom. The main idea is 
that by applying the spline interpolation to the logarithmic and step function terms 
in Eq. (j2- 53|) . we can avoid the situation wherein the discretization points directly 
hit the boundary of the crescent-shape region of the q-q' plane. 

A general prescription is given in §4 of Liu et al.'s paper P^ 1 We first separate the 
q-q' plane into two regions, one is the rectangular region with q < qM and q 1 < qM, 
and the other is the region with q > qM or q' > qM- Here, qM = y <?o 2 — K d 2 (see 
Eq. (|2-52p ). In the latter region, we can prove that xq defined by xq = xo(q,q') = 
((3/4)qM 2 — q 2 — Q ,2 )/(qq') is always less than —1, so that the unsubtracted expression 
in Eq. (|2-53p is safely used. In the rectangular region, the subtraction is made in the 
following scheme. First, let us consider the angular integral 



for q, q' < qM- The kernel function F(q,q';x) is composed of the q, q' depen- 
dence from the spline interpolation for p and p', G(q,q';x), and the rearrangement 
coefficients B(q, q'; x): 




(2-96) 



F(q, q'; x) = G(q, q';x) B(q, q; x) . 



(2-97) 



More explicitly, we move (l/p^ K ) A and (l/p^ K ) A ' factors in Eq. (|2-41|) to G(q, q';x) 
and assign 




The function B(q, q'; x) is expressed as a finite angular-momentum sum 




(2-99) 



k=0 
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with B^ k \q,q') being a polynomial of q and q' . We modify Eq. (|2-96p to 

G(q,q';x) - G(q,q';x ) 



dx 

-l x-x 

, nl / s, f 1 , B(q,q';x) 

+G(q,q;x ) dx - 

J_l X - Xq - lO 



B(q,q';x) 



(2-100) 



Then, the second integral is expressed using the Legendre function of the second 
kind, Qk(z) = (1/2) f\dx Pk{x)/{z — x). Thus, we find 



1 , G{q,q';x) - G(q,q';x ) , 
dx B{q, q ; x) 



i 



X — Xq 



+ (-2) G(q,q';x )Y,(2k + l)B( k \q,q') Q k (x + iO) , (2-101) 



k=0 



where 

Qk(xo + iO) 



Pk(xo) 



log 



1 — Xq 



1 + Xq 



+ iir9(l — \xq\ 



+ 2W k _ 1 (x ] 



We are interested in the integral 







QM 



dq / dq'f(q) I(q,q') g(q') , 
J o 

with the spline interpolation 

/(?) =^2,S IJ ,{q)f{q tl ) for q < q M etc. 



(2-102) 



(2-103) 



(2-104) 



i.e., / = ^2^ u ffj, g u with = ^/U^ f{q ll ) etc. For the integrals to calculate I^ u , 
we can replace the integral variables q and q' with q^ and q v safely if the variation 
of the functions is sufficiently smooth in the mesh-point intervals. Thus, we find 



where 



i rQM rqiu 
—— / dq / dq' S^q) I(q,q') S v {q') 



v^V^ t dx G(^u,x)-G( qil , qv] x 0llv ) B{q ^ x) 

J—l x — Xq^ v 

k M 

Qv] Xq^i/ )^(2fc + l)B (l:) (g„^ Q klw , (2-105) 

fc=0 



Qkfjiu — ( 2) 



1 



U^idv Jq 



QM 



QM 



dq SJq)Q k (,x + i0)S v (q') . 



(2-106) 
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Using these results, we replace the subtracted expression of A^^y^^t iu Eq. (|2-53p 
with 



A 



1 



X()fiv 



PufiK 
Pu^O 



(P^LUK 



B 



~\~Si uu Q S 



(2-107) 



where B^^ v is generated from 5 m7 v, by just modifying Pk{x Qilv ) in g ia >{q^,q u ; 
XQiiv) to Qk^iv in Eq. (|2T06p . Unfortunately, a complete analytical evaluation of Eq. 
(|2- 106j) is not possible. Here, we use the one-side spline interpolation formula and 
give in Appendix B a detailed procedure to calculate 



Qk/iu — ( — 2) 



1 



qu 



dq' Q k (x ^ + iO)S v (q') 



(2-108) 



with xo M = ((3/A)qM 2 — qp 2 —q' 2 )/{q[iq')- Actually, this procedure breaks the symme- 
try of Qkpu with respect to the exchange of /i and v. This generic inaccuracy of the 
spline interpolation technique is, however, very small and we recover the symmetry 
of the S-matrix at the stage of calculating X^s c ^ '5/ in Eq. (|2-62[) . i.e., by modifying 

Xes c ,e'S' c to (l/2)[Xes c ,e'S' c + -^Q's^sJO 

§3. Results and discussion 



3.1. Total cross sections 

For the energies above the deuteron breakup threshold, we further separate the 
momentum intervals over q in Eq. (|2-52p into the following 6 intervals: 



9m 



[0, y M ] 



\qu, 2 ?*f 



^■qu, qu 



[qo,2q ] [2<? ,oo) 



(3-1) 



where qu = y qo 2 — k^ 2 - This is chosen from the criterion that the q^ points do not 
hit the positions of the logarithmic singularities an d th at sufficient points cover the 
rapidly changing region of the kernel in Eq. (|2-107p The first three intervals are 
discretized with the ni-point Gauss-Legendre quadrature, and the next two intervals 
with the n2-point Gauss-Legendre quadrature. For the outermost interval [2*70,00), 
we apply a mapping = 2qo + tan{(7r/4)(l + Xj)} with xi S [—1,1] (i = 1 - n^) 
being the n3-point Gauss-Legendre quadrature. We have altogether 3ni + 2ri2 + n% 
mesh points for the whole q. For the energies below the deuteron breakup threshold, 
we actually use 



<?m : [0, |go] [\qo,qo] teo,2<? ] [2%, 3] [3,6] [6,00) 



(3-2) 



*' Various methods to symmetrize Qk^v cause a serious numerical problem at the point q = 
(V3/2)q M for the solution of Eq. fl^SS}. 

**" ) If we choose the interval [0, qu] with the odd number of Gauss-Legendre quadrature mesh 
points, the middle point hits the logarithmic singularity point (qu/2, qu/2). See Eq. (|B-7|) . 
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where the unit is in fm~ and qo is the incident momentum of the neutron with 
E = (3/i 2 /4M)(/o 2 + £d- Each interval is discretized with the n\-ri2-nj, point Gauss- 
Legendre quadrature in a similar way. Since the maximum value of qo is 0.267 fm -1 , 
2qo does not reach 3 fm -1 . For the bound-state problem of negative energies, we use 

g M : [0,0.2] [0.2,0.5] [0.5,1] [1,3] [3,6] [6,oo) , (3-3) 

although such structure in the small momentum region may not be necessary. For the 
p mesh points, the previous four-interval separation for the bound state problercP 

Pl : [0,1] [1,3] [3,6] [6,oo) , (3-4) 

is used with the ri2-point (for the first three intervals) and ri3-point (for the outer- 
most interval) Gauss-Legendre quadratures. Actual Faddeev calculations, however, 
are carried out using the discretization points only up to 6 fm -1 , to avoid the inaccu- 
racy caused by the spline interpolation. The high-momentum p is, however, necessary 
for the accurate evaluation of the off-shell t-matrix, so that another set of discretiza- 
tion points with 10+15 points is employed to calculate the t-matrix for the positive 
energies. The middle point is the on-shell momentum po = w (3/4)(qa/ 2 — qo 2 ) in 
Eq. (|2-50p . to which the Noyes-Kowalski formalism is again conveniently applied. 
For the negative energies, the discretization points in Eq. (|3-4[) are directly used. 

We also need to consider the discretization points for the angular-momentum 
projection of the rearrangement coefficients in Eqs. (|2-53p and (|2-107p . For energies 
below the breakup threshold, the previous 20-point Gauss-Legendre quadrature for- 
mulsP is safely used for the Legendre polynomials, since there is no singularity point 
for the x integral. For the energies above the breakup threshold, the [—1,1] interval 
for x is separated into two parts, [— l,xo^ u ] and [xo^ u , 1]. For the larger interval, the 
15-point Gauss-Legendre formula is used, and for the smaller interval, the 5-point 
formula is used. 

The three-body model space is mainly specified by the maximum value of the 
two-nucleon angular momentum I max i n Eq. (|2-9p . The maximum orbital angular 
momentum for the two-nucleon subsystem is therefore A max = I max + 1- We need 
to take a sufficient number of the relative angular momentum I corresponding to q, 
to reproduce the backward rise of the differential cross sections. Here, we assume 
4nax = Min{2A max , 10}, which leads to the total angular momentum up to J max = 
4iax + (S c ) max < 3A max - \ since (5 c ) max = I max + \ = A max - \. For the deuteron 
channels, however, the maximum J value is much smaller and J max = £ max + | < 
2A max + § = 2/ max + 1 . We also examine the model space composed of 3 S\ + 3 D± plus 
1 Sq NN interactions only, which corresponds to the so-called five-channel calculation 
of the J 77 = i + bound state. We call this the S + D model space. For example, the 
three-body angular-momentum states in the deuteron channels are restricted up to 
4 Gii , 4 ii5 , 4 Ki9, and 4 Mb, for I max = S + D, 2, 3, and 4, respectively, in the usual 

2 2 2 2 

spectroscopic notation. 

The elastic and breakup total cross sections up to E n = 40 MeV, predicted 
using fss2, are plotted in Fig.[Tl together with the experimental data. Here, the 
incident neutron energy, E n = (3/2)E cra , is measured in the laboratory system. In 
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this calculation, we have used n = ny-ni-n^, = 5-6-5 and / max = 4. Although some 
discrepancy might exist around E n = 10 MeVpD the elastic and breakup total cross 
sections are reasonably reproduced with the constraint of the optical theorem. 

3.2. Differential cross sections 

The differential cross sections for the nd elastic scattering are calculated from 
the scattering amplitudes f^i S ^ ^ s ) in Eq. (|2-92[) by summing up over the final spin 
states and by averaging over the initial spin states: 

S' C ,S C L 

B L (S' C , S c ) = (-lf°- s ° £ Z{W 2 J 2] S' c L)Z{hJ^ 2 J 2 ; S C L) 

x f(e\s> c ),(e 1 s c ) f(e 2 S' c ),(E 2 s c ) ' ( 3 ' 5 ) 

where the rearrangement factor is given by the Wigner-Racah coefficients as^ 1 

Z(£ 1 J 1 e 2 J 2 ; S C L) = l&Jih (£i0£ 2 0|L0) W{h Ji£ 2 J 2 ; S C L) . (3-6) 

The elastic total cross sections are obtained from L = components in Eq. (|3-5|) by 
using a special case B (S' C ,S C ) = X]j«'( 2 ^ + ^)\f(e>s>),{is c )\ 2 - We ^ n< ^ 

= f E B ^ ^ = f E ( 2J + !) E i/<UW a > ( 3 - 7 ) 

S£A (£S C )J (t'SQ 




Fig. 1. fss2 predictions to the nd total cross sections up to E n — 40 MeV, compared with the 
experiment. The experimental data were taken from Ref. I32[| for the filled circles with error 
bars and Refs. I33H36P for the others. 
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which is of course consistent with Eq. ()2-94p . 

We show in Figs. [2] - 0] the nd elastic differential cross sections predicted using 
fss2 for the neutron incident energies from E n = 1 to 65 MeV. Here, we have used 
^max = 4 and n = 5-6-5 to obtain the well-converged values except for E n < 2 MeV. 
(The difference between I max = 2 and / max = 4 is shown in the panel of E n = 3 
MeV as an example.) These are compared with the nd data plotted with bars. The 
experimental pd differential cross sections are also plotted with filled or open circles, 
unless otherwise specified. For the energies E n = 1 and 2 MeV, we find that the nd 
data and pd data are fairly different from each other, while this difference gradually 
diminishes for the energies E n > 3 MeV except for the forward angles 9 cia < 30°. 
We can therefore compare our results with the pd experimental data except for 
this angular region. We find a satisfactory agreement with the experimental data. 
In particular, the agreement around E n = 17 - 22.7 MeV is excellent, which is a 
common feature with the predictions using the meson-exchange potentials!^ For 
higher energies E n > 35 MeV, we find that the forward differential cross sections are 
slightly overestimated in our model. It should be noted that in this energy region, 
many partial waves contribute to Eq. ()3-5p . and yet the shape of the differential cross 
sections is rather simple owing to the strong cancellation. The cross section minima 
(diffraction minima) around 9 cm = 120 - 130° therefore afford a very crucial test 
of the two-nucleon interaction. We find that the minimum values of the differential 
cross sections have an opposite energy dependence to the one given by the standard 
meson-exchange potentials. This energy dependence is very important to discuss the 
effect of the three- nucleon force for the meson-exchange potentials, which is generally 
known as the Sagara discrepancy^* for the disagreement in the diffraction minima 
between experiment and theory. 

In order to investigate the energy dependence of the diffraction minima in more 
detail, we have to incorporate the Coulomb force in our calculation since the precise 
data are only available for the pd scattering. Here, we estimate the Coulomb effect by 
using the published results for the nd and pd cross sections in Ref. 48) for the AV18 
potential. Namely, we use the difference in the nd and pd cross section minima 
in Table II and add it to our calculated results for the nd scattering. The force 
dependence on the difference is considered to be rather small, since the Coulomb 
force is a long-range force. Table [I] shows such a comparison with the experimental 
pd data. We find that on the low-energy side with E n < 5 MeV, our estimated 
values reproduce the experimental data with an accuracy of less than 1 mb, while on 
the high-energy side with E n > 35 MeV, our results are slightly overestimated. As 
for the low-energy side, we will show in a separate paper that the doublet scattering 
length 2 a of the low-energy nd scattering is also consistently reproduced by fss2.'^l | '^3 | 
These results are in accordance with the bound-state calculation of the tritonj^ 1 in 
which fss2 predicts a nearly correct binding energy close to the experiment without 
introducing the three-body force. In the high-energy region with E n > 35 MeV, 
it is reported that the Coulomb effect on the diffraction minima is rather small 
and the inclusion of the A isobar gradually becomes more important to increase 
themPS'E^ These observations imply a possibility that the rather large effect of the 
three-body force, required for all the standard meson-exchange potentials, is related 
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Table I. Comparison of the minimum values of the elastic differential cross sections with the ex- 
perimental pd data after the Coulomb correction. The minimum value (da /d(7) m i n (nd) at the 
minimum point # cm is calculated for the nd scattering at the neutron incident energy E n . The 
Coulomb correction, [(pd) — (nd)]AVi8, for the difference between the nd and pd scattering is 
estimated from the results in Table II of Ref. 48) for the AV18 potential, and "sum" stands for 
our approximate prediction for the pd scattering after the Coulomb correction. The experimen- 
tal data are taken from Refs.l38|) and[48} for E n < 28 MeV, Ref. [46]) for E n = 35 and 46.3 MeV, 
and Ref. S3) for E n = 65 MeV (exp: E n = 64.8 MeV.) 



^ivie v ) 


n 

"cm 

(deg) 


(da J di I ) min(na) 
(mbj 


[(pd) - (nd)jAvi8 
(mbj 


sum 
(mbj 


(da/dn) m y n (pd) 
( m "J 


1 


66 


142.8 


28.6 


171.4 


170.2 ± 1.3 


3 


103 


89.4 


3.3 


92.7 


91.1 ±0.7 


5 


112 


50.3 


3.1 


53.4 


52.7 ±0.4 


7 


116 


31.3 


2.4 


33.7 


32.9 ±0.2 


9 


120 


20.3 


1.8 


22.1 


21.8 ±0.2 


10 


121 


16.5 


1.6 


18.1 


18.0 ±0.2 


12 


123 


11.3 


1.0 


12.3 


12.2 ±0.1 


16 


125 


5.84 


0.5 


6.4 


6.2 ±0.1 


18 


127 


4.47 


0.4 


4.9 


4.7 ± 0.1 


22.7 


128 


2.81 


0.1 


2.9 


2.89 ±0.03 


28 


129 


2.06 





2.1 


2.19 ±0.02 


35 


131 


1.72 






1.52 ±0.04 


46.3 


132 


1.43 






1.28 ±0.02 


65 


133 


1.05 






0.873 ± 0.045 



to the local form of the strong repulsive core, introduced phenomenologically in the 
short-range region. To confirm this, we need to investigate other 3iV observables, 
including the spin polarization and the deuteron breakup processes. We have already 
obtained some good results, especially for the vector-analyzing power of the scattered 
neutronpS'ES w hich we plan to report in a forthcoming paper. 

§4. Summary 

We have applied our quark-model ./V./V interaction fss2 to the neutron-deuteron 
(nd) scattering in the Faddeev formalism for systems of composite particles. The 
energy dependence of the quark-model RGM kernel is eliminated by the standard 
off-shell transformation utilizing the 1/yfN factor, where N is the normalization 
kernel for the two three-quark clusters. This procedure yields an extra nonlocal- 
ity, whose effect is very important to reproduce all the scattering observables below 
En < 65 MeV. In this paper, we have developed our basic formulation to solve the 
Alt-Grassberger-Sandhas (AGS) equations-^' in the momentum representation, us- 
ing the off-shell RGM t-matrix generated from the energy-independent renormalized 
RGM kernel. The Gaussian nonlocal potential constructed from the fss2 is used in 
the isospin basis. 161 The singularity of the iV./V t-matrix from the deuteron pole is 
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handled by the Noyes-Kowalski methodP^'^ Another notorious moving singularity 
of the free three-body Green function is treated by the standard spline interpola- 
tion technique developed by the Bochum-Krakow group Together with the 
results in separate papers^"® discussing the low-energy effective range parame- 
ters and the elastic scattering observables, we have found many new features that 
seem to be related to the characteristic off-shell properties possessed by the quark- 
model baryon-baryon interaction. These include: 1) a large triton binding energy, 
2) reproduction of the doublet scattering length 2 a, 3) energy dependence of the 
diffraction minima of the differential cross sections, and 4) maximum height of the 
nucleon-analyzing power A y {9) in the low-energy region E n < 25 MeV. Further in- 
vestigations on the polarization observables and deuteron breakup processes will be 
discussed in forthcoming papers. 
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Appendix A 

Method to Calculate lC(p,p') 



To calculate W in Eq. (|2-4[) . it is convenient to use the relationship 

W = A(-$=h-^=-h\A = Kh + hlC + lChK, (A-1) 



- ( J=-l).l. (A-2) 



with 



and calculate fC(p,p') in the momentum representation. Since the Born kernel of 
h = ho + Vb + G is already calculated, we can easily obtain W(p,p') by just a simple 
numerical integration. In principle, the kernel K, is calculated from the power series 
expansion in AKA as 



fl2^ {AKA) ' • (A ' 3) 



Quark-Model Baryon-Baryon Interaction Applied to nd Scattering 



25 



where the expansion formula 

1 ^(2r-l)!! 



. 1 = Y 1 /n J" x r for \x\ < 1 , (A-4) 

VT^ fri (2r)H 

is used. In the following, we will show that the infinite sum in Eq. (|A^3j) is taken 
analytically, by using the power-series property for the eigenvalues of the exchange 
normalization kernel K. 

We first consider, for simplicity, a single-channel system with only one quark 
(or nucleon) exchange and at most one Pauli-forbidden state. The Pauli projection 
operator is A = 1 — |uoo)(^oo| with uqq being the h.o. Pauli-forbidden state defined 
in Eq. (|A-12j) below. The normalization kernel in the Bargmann space is expressed 
as 

N = e z *- z ' +X N e TZ '- z ' , (A-5) 

where Xn is the spin-flavor (or spin-isospin) factor and r is given byT = 1 — 1/ /j, with 
\i being the reduced mass number. For the na system, Xn = — 1 and r = —1/4. 
For the (3q)-(3q) system in the QM baryon-baryon interaction, X/y is calculated 
numerically and r = 1/3. We also have to consider the core exchange term in this 
case like in the aa system, which will be discussed later. In the h.o. basis, K is 
expanded as 

K = -X N e™*- z ' = (-X N ) \N)t n (N\ , (A-6) 

N 

where |JV) is the h.o. states with the h.o. quanta N = (N x ,N y ,N z ) and iV = \N\ = 
N x + N y +N z = 2n+l is the principal quantum number. The eigenvalue of K is given 
by 7at = (— Xn)t n (= (—1/4) n for na). The rth power of K is easily calculated as 



K r = (-X N ) r e rrz *- z ' = (-X N ) r ^2\N)(r r ) N (N\ 

N 

= (-X N ) r \N£m){T r ) N '(N£m\ . (A-7) 

Nlm 

The kernel for [AKA) r can be obtained by restricting the sum in Eq. (|A-7j) over 
N> 1. 

Suppose the spatial part of the GCM kernel for K is i N (z-,z') = e TZ "- z '. The 
corresponding Born kernel is given by 

2vr 1 \ 3/2 f 1 (1-t 2 1 + t1 2 



J2(2l+l)M e (q f ,q l )P e (q r q l ) 

47T^M,(g / ,%)^ m (q / )y/ m (q,) , (A-8) 



(in 
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where k = qf — q i} q = (q^ + Qj)/2, and 7 = \xv with v being the h.o. width 
parameter of clusters. (See Appendix A of Ref. l29|) ,) Thus, for the partial wave 
component 2f| pace — Y^^Lo \N£m)T N (N£m\ with N = 2n + £, we find 



3/2 



2tt 1 
7 1-r 5 



3/2 



exp 



1 + 1 ( 2 , /2 , 



1 



1 - T 2 47 

3/2 ,2 

e 



1 r 



7 1 — T 2 



T- 8 1 



(, 2 + g ' 2 ) /l r 



1 9 

7 1 — r z 

where is the imaginary spherical Bessel function. We use the notation 

x 



(A-9) 



i e (x) = i e ji(-ix) 



(2£ + l)W 

eK ) ^ (2n)!!(2n + 2^ + 1) 
1 



r-FHx 2 



ra=0 
1 + 



-X + 



1 



2- (2^ + 3) 2-4 - (2£ + 3)(2^ + 5) 
By using this notation, Eq. (|A-9j) can be expressed as 

1 \ £+3/2 



X 2 + 



(A40) 



Me(q,q']T) = t u £{q)u oi (q') 



1-r 2 
/\ 2 \ 



1 — r 2 7 



where 



wo«(g) = u e(q;j) 



2vr\ 



3/4 



1 

e 4 t 



(A-ll) 



(A-12) 



7 / V( 2£ + 1 ) !! 

is the h.o. wave function with I and the lowest h.o. quanta n = 0, normalized as 



47T 



^3 y g 2 <ig (uoi(q)) 2 = 1 . (A-13) 
Thus, we find 

00 

K\ = {-X N f (r r ) 2n+e u ne (q) Une (q') = (-X N yM e (q,q';T r ) . (A-14) 



n=0 



Here, we treat only the n = term separately. Namely, by defining a new function 



Mi(q,q';T) 
we find 



1 \ 



£+3/2 



1-r 2 



_T_qq_ 

1 — r 2 7 



tfj = (-A^r^^H^') 1 + M £ (g, </; r r ) 



1 , (A-15) 



(A46) 
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Thus, if there exists no Pauli-forbidden state, A = 1 and 



r=l 



(2r- 1)!! 
(2rV 



-Kl 



Or — 1 \\\ 

E ^J" (-X N rrnodq)u e(q') 1 + M t (q, q'; r r ) 



r=l 



(2r)H 



(A47) 



For the first term, we take the r sum with Eq. (|A-4p . Then, we eventually obtain 
JC e (q,q') = u oe (q)u oe (q') 

(2r-l)!! 



V 1 + x N r e 



\ (2r)! 



:(-X^rM^(g,g';r r ) 



(A-18) 



If a Pauli-forbidden state exists only for (Os) with £ = (namely, 1 + Xn = 0), we 
have Xnt £ = — 1 and the first term of Eq. (|A-17j) should be omitted. Namely, 



£o{q,q') = E 



" (2rJ\ ( AKA >t=o 



™ (Or - IV 

M<?)*W) E (2r)!! " Mo(g ' ^ ^ 

r=l ^ 



(A-19) 



For £ 7^ 0, the subtraction of one in Mi(q, q'; r) seems to be redundant since the 
(— 1) term cancels with the first term in Eq. (|A-16j) . However, the convergence of 



v/i + x NT * 



i = E 



(2r- 1)!! 
i (2r)!! 



(A-20) 



is very slow if (— XjqT^) is close to 1. This happens in the P-wave case of the QM 
NN interaction. Namely, for the NN interaction, we find (see Table II of Ref. l5ip ) 



(A-21) 



On the other hand, M^q, q'; r) is a function of r 2 with Mp(q, q'; 0) = 0. The leading 
term in the power series expansion of r 2 is 



state 


Xn 




IN 


-V 


Si, 1 5o 


i 

9 




1 (l\ 2n 
9 \3J 


1 
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'Pi 


7 
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/l\2n+l 
1 37 


7 
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3 Pi 


31 
27 


31 
27 




31 
81 



M £ (q,q';r)=T 2 



+ l-^(q 2 + q' 2 )+ ' 



2 2 7 



2(2^ + 3) V 7 



/\ 2 



+ 



(A-22) 



Thus, we have M £ (q,q';T) = 0(r 2 ) and M e (q,q';T r ) = 0(r 2r ). Since r 2 = 1/9 in 
the above example, increasing r gives very fast convergence on the order of 1/10 for 
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r = 1, 1/100 for r = 2, • • • . We take the r values up to rMax = 15 in the actual 
calculation. 

In the application to the (3q)-(3q) RGM kernel, we need a proper treatment of 
the core exchange term and the coupled-channel problem. In the operator formalism 
of the spin-flavor factors, the basic (3q)-(3q) GCM normalization kernel is expressed 
as 

I%(z; z') = e z *- z ' + X N eh z *- z ' , (A-23) 
by which the full normalization kernel is given by 

N{z-z') = \ ll%{z-z')-P a P F I%{z--z')} . (A-24) 

Here, Pp is the flavor exchange operator. Thus, the exchange normalization kernel 
K is expressed as 

K(z; z') = ~ {X N e\ z "- z ' - P a P F X N e~* z *- z ') . (A-25) 

The spin-flavor-color factor Xn contains an explicit Pp dependence, if the bra and 
ket sides are composed of nonidentical baryons: 

X N = X d + X e P a P F . (A-26) 

The exchange operator P a Pp takes the value —(—1)^, where (—1)^ is the parity of 
the two-baryon system. Thus, Xn has an explicit parity dependence 

X £ N = X d - {-lfX e , (A-27) 

which is important in actual calculations. The partial-wave component Kg is then 
given by 

Kg = -Xi r J2\N£)U) (Nt\- (A-28) 

From now on, we omit the superscript i in X l N , by assuming a fixed I. The eigenvalue 
problem of the multichannel Kg is reduced to the eigenvalue problem of the matrix 
(X N ) a p. We solve 

T,( X »U C fi = XC Z> (A-29) 

P 

with 

a X 

Then, the exchange norm kernel is given by 

N 



(Kg) a p = -Y,^C*\N£)(\) (N£\. (A-31) 



Quark-Model Baryon-Baryon Interaction Applied to nd Scattering 29 

The full eigenvalue is 7an = — A(l/3) Ar . Only the (Os) state is possible for the 
Pauli-forbidden state in the isospin basis; i.e., the SU$ (ll) s state: 



{AK,A) ap = £ (-A) Clcl \Nl) 



N 



(N£\ 



((AK £ A) r ) ap = (- A ) r (3 J 



<iW| 



^(^) a/3 -«5, j0 CiCi|00}{00| . 



(A-32) 



AiV 



The rest is almost the same as in the single-channel case. The final result is 



(2r-l)!! 



E 



— 1 I C^Ca 



+E 



- (2r)!! 



?(-«)') 



3 



(A-33) 



Appendix B 

Method to Calculate Qk^u i> n Eq. \2-108\) 



We first separate the integral region of Eq. (|2-1U8|) as 

Qkjxv = (-2) — V] / dq Qk(xo fl + iO)S u (q) , 
with qo = and q KM = qu, and apply the third-order spline function 



(B-l) 



#(5)=E^ H (?-««r for ge[? K -i,« K ]. (B-2) 



m=0 



Then, we find 



K M 3 



-1 ' °m 

Qk»» = (-2)- 52 Yl < M > (b-3) 



k=1 m=0 



with 



Qffi(9«-l,?»e) 



dg' Qk{x 0fl + iO)(q' - q K ) r 



Qk-1 



(B-4) 
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For k = 0, a completely analytical calculation is possible by using the integral formula 
for logarithmic functions 



/ 



dx log(x + q) = (x + q) [log(x + q) — 1] 



dx x m logx 



1 



„m+l 



m + l 



logx 



1 



m + 1 



(B-5) 



and 

I m (a,b; q) 
+ [(b-a 

m—l 



.m+l 



d<?' log 

q + q 

■{b-q) m+1 ] lo. 



(</-&) m = (-l) 
q — a 



1 



(6-g) m+1 log 



m + 1 

[(6 + g) m+1 - (6-g) m+1 ] log 



b-q 



r=0 



b + q 
q + a 

q + b, 

(B-6) 



with m = 0, 1, 2, 



• and < a < 6. We can prove that q > does not need to 
be < a < q < b, but can also be at any position in the above form. The formula 
in Eq. (|B-6p is free from the logarithmic singularities since lim e _>o e loge = 0. We 
define (^-dependent variables, q±^ and q 2fM , for the crescent-shape region: 



qifi 



\Jl(qM 2 -q^ 2 ) - \q^ 
Q2» = \J\iqM 2 - q^ 2 ) + \q^ 



qi» = \q» - \j\(qM 2 - q 2 ) 
q2fM = \q» + J\{qM 2 - g M 2 ) 



for q^ < -^-QM , 



, y/3 

for ~2~qM < Qfj, < qu 



(B-7) 



Then, we obtain 



+iTrI^ l (q K - 1 ,q K ;q lfl ,q 2fi )] for < ^ - }t 

{ ^qu <q»< qu 



(B- 



where the theta function part is given by 

1 



I m (a,b; qi,q 2 



(Mm{b,q 2 }-b) 



m+l 



(Max{a,gi} - b) 



m+l 



m + l L 

for Min{6, q 2 } >Max{a,qi} , otherwise 



(B-9) 



When k > 1, various methods are used to calculate Qmji(qn-i, q K ) in Eq. (|B-4[) 
accurately. First, the most accurate calculation outside the crescent area and its 
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neighborhood is maybe the numerical integration using the power series expansion 
of Qk{x) in 1/x 2 : 

2 k ^ ( n + k)\(2n + k)\ ( 1 \ n c ,„ . 

Q fc ( s ) = — - V ^-77 ^ rr [ -o for X > 1 . (B-10) 

n=0 

The convergence is so rapid that we can use Eq. (|B-10p for the numerical integration 

of Qmi(qK-i, q K ) with x 0fJ ,(q' = q K ) > 1-05 or x 0fl (q' = q K -\) < -1.05. We use the 20- 
point Gauss-Legendre quadrature for the numerical integration. (Note that there is 
no imaginary part appearing in this case.) In the crescent area and its neighborhood, 
we will expand Pk(xo^) and Wk-i(xo^) in Eq. (|2-102p around q' = q K : 

m M 



m'=0 
m M 

^iW=EC M («'-^ m '. (B-ii) 

m'=0 

for </ S [^k-i, <7 K ]. Then, we can use the formulas in Eqs. ()B-6p and (|B-9j) to obtain 

(1 \ mM i 
' m'=0 

+iiri^ +m r(q K -i,qK;qi^q2ix)] +2/?^ K(m) C+ m '(g#t-i, <? K ; <?«)} , (B42) 



for < (\/3/2)g^ and (v3/2)<Zm < q^ < qM, respectively. The expansion coeffi- 
cients /3^ m ^ etc. are expressed using Bell's polynomials!^ 

Y = f , 

Yi = fi9i, 

Y% = fi 92 + h 9i 2 , 

Y 3 = h93 + h (3g 2 5i) + / 3 5i 3 , 

Y± = fi94 + h (%3 91 + 3 52 2 ) + h (652 9i 2 ) + h 9x , • • • , (B-13) 

where the subscripts imply the higher derivatives. If we write Eq. ()B-13j) as Y m (f, g), 
^ku^ e ^ c - are expressed as 

fik\i ^ = ~\ Ym (Pk(x0fi K ), 9^k) , 

C K(m) = ^ Y ™ ^iW.S^) • (B-14) 
The higher derivatives of g = g^ are given by 



<W ' q K ° ' ' ' q K ' 



3i = -(7 i 5 H ) ' 92 = x lx —£, ••• , c/ r = a^(-l) r — , (B-15) 
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with Xfj, = ((3/A)qM 2 — Q^ 2 )/q^- For the practical calculation, we first expand P k (xon) 
etc. around the middle point (q K -i + (/ K )/2 and then rearrange it to the form of Eq. 
(|B-lip . by using itim = 3. Actually, Eq. (|B-12p cannot be used if q K is small. This is 
because the higher derivative of g in Eq. (|B-15j) becomes very large for the small q K . 
This method is not valid either when q^ is small and xq^ rapidly changes from — 1 
to 1. We therefore restrict the use of this method to the region with q^ > 0.6 fm -1 
and q K > 0.6 fm -1 . 

The third method to cover the above missing area is to use the separation of 
Qk(x + iO) to the singular and nonsingular parts: 

Q k (x + iO) = Q k (x + iO) + P k (j^\ Q (x + i0). (B-16) 

Note that Qo(x + i0) =0 and the nonsingular function Q k (x + iO) satisfies the same 
symmetry relation as Q k (x + iO), i.e., Q k {— x + iO) = (—l) k+1 Q k (x + iO) for real x. 
We calculate 

flK _ 

Q%l(q K -l,q K )= / dq' Q k {x 0ll +i0){q' -q K ) m 
J q K -i 

+ J' K dq' Pk(j^j Qv{x^ + iO){q'-q K ) m (B-17) 

separately. The numerical integration is used for the first integral since the integrand 
is nonsingular. However, if the interval [q K -i,q K ] contains qi^ or q^, we separate 
the integral region into two parts, [q K —i, qi^] and [qif_i,q K ], etc. The second integral 
in Eq. (|BT7p is reduced to the previous formula, resulting in 



j dq P k ^j^i^ Qo{^On + i0)(q' - q K 



in 



\\xon\J 

) for x 0f , K > 

(-l) fc Qmi{q K -i,q K ) for atytfc-l < 

Qm]i(9«-1,9«) + [(-I)*"' - 1] Qm|t(?0Ai)?«) for X p K < < X 0flK - 



(B-18) 



for q^ < (\/3/2)qM, where Qm^{a,b) is obtained from Eq. (jB-8[) by replacing q K -i 
with a and q K with b. In the last term in Eq. (|B-18j) . we have defined qo^ = 
(3/4)qM 2 — q^ 2 - In the case of (y/3/2)qM < q^ < qM, only the second case of 
Eq. ([B^T8j) is realized. 
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